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ON THE GENERAL PERTURBATIONS OF THE 
POSITION VECTORS OF A PLANETARY SYSTEM 


by 

Peter Musen 

Goddard Space Flight Center 
SUMMARY 

A theory of the general perturbations of a planetary system is developed in this 
article. 

The perturbations of the position vector of each planet are expanded into series 
arranged in powers and products of the masses mi , m 2 , * • ■ , m n of the planets consti- 
tuting the system. Perturbations of different orders are obtained in the form of series 
containing the purely periodic, the secular, and the mixed terms in accordance with 
standard astronomical practice. The influence of the lower order perturbations on the 
other ones is determined. 

Typical differential equations are formed to determine those perturbations of the 
i th planet which are proportional to m k , m k m p , m k m p m s , ••• . The right sides of these 
differential equations are obtained as the corresponding terms in the Maxwellian expan- 
sion of the gravitational forces in terms of multipoles. The momenta of these multipoles 
are the perturbations of all possible orders. 

The explicit calculation is carried out here for the perturbations of the first, second, 
and third orders; and the procedure for determining the higher order perturbations is 
outlined. 

Decomposing each perturbation of any particular planet m. along the undisturbed 
position vector r. , along the undisturbed velocity v., and along the unit vector R. nor- 
mal to the undisturbed orbital plane, we reduce the differential equations to a form easily 
integrable by quadratures. After the integration, it is more convenient in practical appli- 
cations to replace this decomposition of perturbations by the decomposition along r., 
R^r. , and r. . The problem of the constants of integration is treated for the case of the 
mean elements. The results given here extend and generalize the author* s previous re- 
sults to the case of the whole planetary system. The method suffers, however, from dis- 
advantages common to all astronomical methods of the general planetary perturbations: 
It is not applicable to a pair of planets if their orbits approach each other very closely. 
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ON THE GENERAL PERTURBATIONS OF THE 
POSITION VECTORS OF A PLANETARY SYSTEM* 


by 

Peter Musen 

Goddard Space Flight Center 


INTRODUCTION 

A theory of the general perturbations of a planetary system is developed. The perturbations 
in the position vector of each planet are developed into a series in powers and products of the 
disturbing masses and into a series containing the periodic, the secular, and the mixed terms with 
respect to time. 

Such a way of representing the integrals of the disturbed motion is in accordance with stand- 
ard astronomical practice. From the purely mathematical standpoint, this solution can be affected 
by all the difficulties associated with the near -resonance conditions caused by the small divisors. 

We establish the differential equations for perturbations proportional to the powers and prod- 
ucts of masses in a form integrable by quadratures. The explicit calculation is carried out through 
the perturbations of the third order. In our planetary system it is rarely necessary to include the 
perturbations of the fourth and higher orders. However, an outline of the procedure for including 
the perturbations of these higher orders is indicated here. 

The problem of a direct determining of the general perturbations in the position vectors, 
including the effects of higher orders, became possible only in recent years with the advent of 
electronic computers. By decomposing the perturbations S?. along the directions of r., v. , and 
R. (Reference 1), we can integrate the variational equation of the problem by Hill T s (Reference 2) 
procedure directly without resorting to the method of variation of astronomical constants. We 
shall use here the same decomposition as an intermediary step; but the final decomposition of the 
perturbations will be along r. , R.xf , and R., to reduce the components of the disturbing term on 
the right side of the variational equation to a simple form. In computing the higher order pertur- 
bations, it will be necessary to expand the disturbing forces in powers of perturbations in the 
position vectors. Maxwell’s method of expanding the electrostatic potential in terms of multipoles 
by employing symbolic operators (Reference 3) can be used profitably also in planetary theories. 

In our exposition the moments of the multipoles are the perturbations of different orders of the 
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position vectors. Evidently any other way of expanding the disturbing forces in terms of the per- 
turbations of the position vectors will lead to a duplication of Maxwell's expansion, but through a 
more laborious writing. 

In the theory of perturbations of the position vectors, economy of theoretical thinking as well 
as economy of computing machine time is achieved because a set of homogeneous operations is 
being constantly repeated. All these circumstances suggest that future methods of calculating 
general perturbations will be based on the expansion of the perturbations directly in the position 
vectors. 

THE DIFFERENTIAL EQUATIONS OF THE PROBLEM 

Putting 


p = r - a 


we shall make use of the Maxwellian expansion of the spherical functions as defined in terms of 
multipoles. We have 


<£("> = 




P 


= (-l) n 


Z>-.TTv* 


TT-v?- 


(k= 1,2, ••• ,n; n= 1,2,3, •••)■ 


( 1 ) 


where a,, a 2 , • • • , a n are constant vectors, V is the del operator with respect to 7, and the sums 



IT v?> etc - 

k = 3 


designate the sums of all terms as obtained from the first term by means of the permutations of 
all n indices. In particular, we have 



( 2 ) 


2 


(3) 


1 a , ' P 

0( 1 ) = a • V i = L 

^ 1 o 


p P 3 


0(2) = V V * 2 -vl = t^p-a, P'S 2 --j V*a 

V pS p 3 


(4) 


0(3) = a -Va 2 -Va 3 -V ^ ~ P • a, P • 3 2 P * a 3 + .1 (P • a, a 2 • a 3 + p-2 2 + p • a 3 a,-^), ( 5 ) 

F P 7 P 5 

0( 4 > = 3j -V a 2 • V a 3 • V a 4 - V i = + p*a x p ■ a 2 p • a 3 p • a 4 

P p 9 

- 1£ (p • a 2 p • a 3 a t • a 4 + p • a 3 p • a t a 2 • a 4 + p • a x p • a 2 a 3 • a 4 
P 7 


P • a, P • a 4 3 2 • a 3 + p • S 2 p • a 4 • S 3 + p • S 3 p • a„ a x • a 2 ) 


+ ( a j ‘ a 4 a 2 -a 3 + a 2 -a 4 a^ + a 3 


•3 4 S,-a 2 ). 


(6) 


The gradient of the spherical function 0 (n > is obtained from the Maxwellian expansion of 
0 < n+1 > simply by replacing the moment a n+1 by the idemfactor I. Thus, from Equations 2 to 6 , 


V0(°)= vl = -JL, 

p p 3 


(7) 


- „ _1 3 — ~ - 1 

a r vv p = +-pp- a , - — v 

^ p 3 p 3 


(8) 


V V a 2 -V vl = - Hp p-5 t p-a 2 + -1 (p • *, a 2 + p • a 2 a, H-pa^a^, 


(9) 


J x -V a 2 -V a 3 -V vl = + l^p p • a, p • S 2 p-a 3 
P p7 


15 


(P * a 2 p-a 3 a 3 + p • a 3 p • a t a 2 + p • a x p • a 2 a 3 + p p • a x a 2 -a 3 


+ p p • a 2 a x -a 3 + p p-a 3 a,-^) + JL (a x V S 3 + a 2 V 3 i + a 3 a, -a 2 ) • (1()) 
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The differential equation of the disturbed motion of the i th 


planet can be written in the form 


^ ( r . + Sr. ) = V. 

dt’ 1 11 




S r ; I 


L 

cr~ 1 


fm V (_ 


lp t ri+ S ^J 


r + 


Sr | 

cr 1 


(i, cr - 1, 2, n). (11) 


Taking the equation 



into account and introducing the vectorial differential operators 


( 12 ) 

(13) 

which perform the development of Equation 11 into a power series in the components of 
$r k (k = l, 2, • • • ,n), we obtain from 11 the differential equation for Sr . in the form 


D. = V. exp ( 8 7 . • V. ) , 
D. . = V. exp(S^. . * V. ) 



(14) 


The perturbation vector Sr . can be developed into a series with respect to the powers and 
products of the disturbing masses. We put 


Sr. 

i 



a 


o./3 a/iry 


(15) 


where 7“ is proportional to m a and r, a/3 is proportional to m a etc. The factors in front of 
the sums in Equation 15 are introduced to remove large coefficients in higher approximations. 

We define 7“^, 7. ai37 , • • • in such a way that they remain invariant under the permutations of 
the upper indices: 

7?£ = 7?“ , = 7^ a = ... etc. 

Ill » 

Symbolically, 
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Sr. = f , exp 2^ 7“j - 


1 . 


where in performing the development and the symbolic "multiplications" the indices are not being 
added, but written in a row. We also have 


l _ r lll _ 


= o. 


We now shall deduce the differential equations for determining perturbations of the form 


7 k r kp r kpq 7 kpqs 
i ’ i * i ’ i * 


(i, k, p, • • • = 1, 2, • * * , n) , 


first under the assumption that there are no identical indices among k, p, q, • • • . Retaining only 
the substantial terms, we have 

87. = (7 k + ?P + 7? + 7») + (7** + T 1 ?* 1 + 7 ks + rf + 7 p s + 7 qs ) + (7^ + 7 k P* + 7f* s + 7 m*) 

+ r*M« + ••• , ( 16 ) 

• 8 ?.. = ( ^. + pf L +pf i ) + (p'‘ ? + + + + ••• • (i 7 ) 

Substituting these values into 

87. -V. (87. -V.) 2 

D. = V. + — 1 1 V. + 1 -V, +•••, 


1 1 


1 ! 


2! 


S3.. -7. (8p..-V.)2 

D =7. + 11 J V. + 31 J - V. + ••• 

H j + 1! J 2! J 

and again retaining only the necessary terms, we deduce: 

D. - V. = 7 k *V. V. + (r kp - V. V. + 7 k *V. 7 P *V. V.) 

i i iii v i ii'iiiii 7 

+ (7 l ?P < J-V. V. + 7 k P- V. 7?-V. V. + rP^-V. 7 k -V. V. + rf-7 7 p -V V + 7 k -V rf -V 7 q -V V ) 

v i ii T i i i iii 1111 1 1111 1111111 ' 

F-kpqs.y V + r kpq *V r s -V V + r kps *V. r q -V. V. + 7 kqs .y ?P.y y + r pqs -V. 7 k *V. V. 

' i ii i iiii i 1111 l iiii 1 1111 

/^kp.y ^qs.y y. + ykq.y y P s.y y fks.y ^pq .y V.)+(r k P-V. r?*7. r s -V. V. 

i l ii^i v i 1 i ii^i i i ii /v i 111111 

+ r kq *V 7 s -WFp-y y + r ks *V. r p -V. r q -V. V. + r pq *V. r k *V. rf-7. V. 

r i i i i i i i i i i i a l 1 1 1 i 1 l 1 

7PS.V. 7 k -V. 7 q -V. V. + 7 q . s -V 7 k -v 7P-V. V ) + 7 k -v. 7?-V. 7 q -V. 7?-V. v] + ... fi q\ 

i 1 i 1 1 1 1 1 v ii i i i i * * * 1 1 1 1 i i ij ^ > (lo) 
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( 19 ) 


D.. = V. + p k .V. V. + (p k p-V. V. + p k .-V. p?.-V. V.) 

i> j i j j r ii j r ji i j' 

+ (p k ? q -V. V. + p k P-V. p'l.-V. V. + pP q -V. p k .-V. V. 

J I r l‘ J 'j » J J T ^Ji J 'l i j j 

+ p qk -V. pP.-V. 7. + p k . -V. pP. -V. p q .-V. V.) + ••• . 

These developments of the operators D j - v . and D. . permit computation of the general plane- 
tary perturbations up to the fourth order if necessary. 

We will establish the differential equations for determining the perturbations up to the third 
order. In our solar system occasions requiring the perturbations of the fourth order probably 
will be very rare. However, in some cases of very sharp commensur abilities of mean motions, 
the question remains open and further numerical investigations are necessary. Substituting 
Equations 18 and 19 into 14 and retaining only the typical operators, we have 

— ( r k + r kp + r k P<* + ...)= |~r k -V. V. + (r k P-V. V. + r k *V. 7 p*V. V.) + -V. V. + r kp *V. V. 

j2 V I l l /r xL.XXX TV X X 1 1 1 X 1 X X 1 X XX XX 

dt 


+ r^* 


pq. V ?k.y y 7 <jk.y 7 ? .y y -k.y -p.y -q.y y N .] J_ 

X 1 X 1 X T X xxxx^xxxxxxx 7 Jr. 


+ K {- [v k + ^'V k v k + + •••] ~ 

+ K + 7 r^ v k + ( 7 r- v k v k v k> + •••] 

+ fm p {- + •••] — 

+ rr k -V V + fr qk -V V + 7 q -V r k -7 V ) + . . .1 JL l 

j_PPP V P PP PPPPP' J r p j 


fm 


- [(pkp-V V + p k . -V pP. -V V ) + ...1 _J_ 
[} r q X q Q 'q * q q q / J n 

'q x 

T/^kp.y y 7 k. y ?p.y V ) + . . -1 — L . 
L q qq qqqqq 7 J r f 


( 20 ) 


Comparing the terms of the same degree in the disturbing masses in the left and right sides of 
Equation 20, we obtain the basic equations for determining the general perturbations up to the 
third order: 


d 2 7 k 
1 

dt 2 



v. 

X 



(21) 
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d 2 7 k P 


Fi = /x? fr k P-V. V. L + F kp \ 

dt 2 1 \ 1 ‘ ‘ r i / 


(22) 


d 2 ? kpcl 


_ll _= M 2 (r kpq -V.V 1 + Fi"V 
dt 2 V i / 


(23) 


where 


F k = k V * + V -A-\ , 

1 1 + ra i \ k ^i k r J 


F kp = + 7 k -V. ?P-V. V. — 

* T ii iixr. 


(24) 


™ k /l pP -V V _L + 7 P -V V -L.\ + mp p k . -V V — + 7 k -V V 1 
1 +m, \ Pki k k p,. ; k k k r k j + 1 + m. V Ppl p p P„, p p p r 


(25) 


and 


F k P q = (r k P-V 7 q -V V + rP q -V. 7 k -V. V. + 7 <)k -V. 7 P -V. V. + 7 k -V. 7 P -V. 7 q -V. V. ) — 

r i v i iiiii l l l l T i liiiiiliiii'r. 


1 + m. 


n<i k 


1 + m. 


- (p qk -v v + p q .-v P k .-v vj _L + (7 qk -v v + 7 q -v 7 k -v v ) J_ 

vr pi p p pi p P 1 PP ,o. v p pp p p p p p 7 r 


1 + m. 


- (p k P-V V + p k . -V pP. -V V ) _L + (7 kp -V V + 7 k -V 7 p-V V ) -L 

' ' q i q q r qi q r qi q q 7 n v q q q q q q q q r 

' q i 


(26) 


Equations 21 to 23 can be written in the form 


d 2 7 k / i 37. 7 A 




•7 k = m 2 F k 

i ~ l l 


(27) 


d 2 7 kp / T 3 r. r \ „ 

i_ +M 2 (— 4— )' r i P = Pi Fi 

dt 2 ' ’ ' ' 


kp 


r? r? 

V i 1 


(28) 


d 2 ? kpq / 1 37. 7. A 

+ M i(— - — : — ) - 7 * 5^0 = /i 2 F- pq - 

l ~ l 1 


dt 2 


r? r 5 


(29) 


The terms in the right sides of Equations 27 to 29 are the partial gradients of the sums composed 
of the elementary spherical functions, with the moments equal to the perturbations in r. and 
p (i, k = 1, 2, ■ , n) . Making use of Equations 24 to 26, we obtain the following: 
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(30) 


F k =- 


Pn r k 


1 + m s \ o 3 . 


k i k / 


also the expanded forms of the typical terms in F. kp : 


7 k *V. 7?*V. V. ±. = 

i i x i i r . 


15 


. r . *7 k r. *r p + (r . *F k r p + r. *F P F k + F. r k ‘r p ), 

ii ii i , v i i ii i ii i i /s 


•F k F p 


(31) 


mk p P -V. V. 1 


1 + m. ' ki 1 1 p k . 1 + m i l ^3 


k I _I ^ i Pki 

Pki 


Pki 


(32) 


mk _ „ 1 

+ -r Fj-V.V. — = 

1 + m. k k k r 


n k I 


3 r k r k 


1 + m i Vr? r? 


(33) 


and the expanded forms of the typical terms in 


• kp -V. 7 q ■ V. V. — 7. 7. -7 kp 7. -7 q + (7. *7 kp 7 q + r . •r < ? ? kp + 7 7 kp, ? q ' 

i l 1 1 1 r . n ill 11 C V 11 l'l i i T i i i y 

1 r 7 r? 


(34) 


? k -t7 7 p -v 7» -V. V. — = + i££ 7. 7. -7 k 7. -7 p 7. -7 q 

-ii li ill 


i i i i i i 1 r . 


_ ^ (7. ■?? r •?*? 7 k + 7. *7? 7. -7 k 7 P + 7. *7 k 7. *7? 7? + 7. 7. *7 k 

— - \ i i i i i i i i i i i i i i i T i i i 

+ 7. 7. -7? 7 k -7*? +7. 7. ■?<? 7 k *7?) 

ii i i iii i i i / 


+ _L (7 k 7 p -7 < ! + 7 p 7 q -7 k + 7 q 7 k -7 p ), 

c v i i i i i l l l l ' ’ 

rr 


n i 1 

k 7- = ■ 


1 + m . ■ kl K K Pki 1 + m i\pki 


3 Pki Pki 


Pki 


m, 


1 + m. 


p p. -v. pi. -v. v. _J_ = _IL 

k ‘ k ^ kl k k / 3 k . 1 +m . 


— — Pki P\l i "P k i Pk i "P k i 

Pki 


+ — (Pki'PkiPki +Pki - Pki Pki +Pki Pk P i-Pki) 

Pki 


(35) 

(36) 


(37) 


n. 1 

k 7 pq -V V — = 


I 3 r k V 


1 + m. k k k r k 1 +m. I r 3 


k ’ 


(38) 


rn. i 

— — 7? -v. 7? *v u v — - = 

1 + m . k k k k kr k 


1 + m. 


15 


r k 7 k’ r k V r k -■ 


(?.•?£ 7? + 7,-7; 7 p + 7 

v k k k k k k 


. p -7 q ) 

k k r k' 


(39) 
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In a similar way the expressions for F^ pqs and for its typical terms can be formed. The process 
can be continued as far as necessary, until the perturbations become negligible. The restriction 
that indices k, p, q, be different now can be removed. 


INTEGRATION PROCEDURE 

The typical differential equation which appears in the problem of the general perturbations 
in the position vectors has the form 


d 2 X. /x. 3 X •?. r \ 


(40) 


where 


tp _ pk ijp k P ^ k pq .... 

r i “ r i ’ 1 i ’ ' i * 

Fi is a series with the periodic, secular, and mixed terms. Equation 40 is the variational 
equation of the two-body problems. 

To reduce the solutions of Equation 40 to quadratures, we shall make use of the substitution 

x. = (s. + 2w.) r . - v. J (2s. + 3 w. ) dt + £. R. . ( 41 ) 

This differs from the substitution used by the author (Reference 1) previously. The substitution 
of Equation 41 was chosen because of a simpler form to which Equation 40 is reduced thereby as 
compared with the earlier exposition. It is unnecessary to retain the index i in the further 
exposition; we can now omit it without loss of clarity. 

It follows from Equation 41 that 

§=fe* 2 t ,42) 


d 2 x 

d 2 s 

0 d 2 w 3/j 2 
+ 2 - ___ 

dr 

(2s + 3w) dt + t— (3s + 4w) r + 

£ + ft f (2s + 3w) dt 

dt 2 

dt 2 

dt 2 r 4 

dt 

1 r 3 J 

L dt r 3 J J 


v + £!£ R- (43) 

dt 2 


Substituting expressions 41 and 43 into 40, we obtain the vectorial differential equation 


d 2 s d 2 w\ -* dw — / d 2 £ (m 2 r \ — 2^ 

+ — s + 2 I r + — v + + J — c, R = M F , 

3 dt 2 / dt \ dt a *» 3 


dt 2 


(44) 
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which can be integrated by quadratures. Forming the dot products of Equation 44 and 

vxR, Rxlr, R 


and taking 


r vxR - v*Rx 7 - M ]/a ( 1 - e 2 ) , 


r 


d r 
dt 


into account, we have 


d 2 s 

d7 2 



s + 2 


d 2 w 

dt 2 


na 

YT^ 


F • vx R , 


(45) 


dw 

~di 


- na — F-Rxr , 

YT^ 


(46) 


dH, 

d t 2 


+ — C = /u 2 F-R ■ 


We have from Equation 46 


where 


w = k 3 + b , 


B = 


I 


na 

/l - e 2 


F*Rxr dt ; 


(47) 


(48) 


(49) 


k 3 is the additive constant of integration. The integration of series in Equation 49 is performed 
in a formal manner. 

Equation 45 now can be integrated by using Hill's procedure (Reference 2). We obtain 


s = K !~ cos f + K — sin f + A, 
<a * a 


(50) 


where K x and K 2 are constants of integration and 


= 1 


(f-vxr) (R.?xf) 
a (l-e 2 ) 


dt - 


I' 


J2 d 2 w 

dt 2 


( R-rx r ) d t , 


(51) 


The vector f is considered as a temporary constant and is replaced by r after the integration is 
completed. The integrand is a trigonometrical series in the mean anomalies of planets and the 
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auxiliary mean anomaly l associated with 7 ; it also can contain the purely secular and the mixed 
terms. After the integration is performed, 1 has to be replaced by l . Integrating the second 
integral in Equation 51 by parts and replacing r*7 by zero when it appears outside the integral 
sign, we obtain 


f^R.?x?dt = -f^R.vx 
J dt 2 J dt 


7 dt . 


From Equations 46, 51, and 52 we deduce 

A = f : (f.vxR) dt + f (f-Rx? ) 

J a ( 1 - e 2 ) J a (1-e 2 


dt . 


As in Equation 49, the integration is performed in a formal manner. 
In a similar way we obtain from Equation 47 


(52) 


(53) 


l = K c — cos f + K fi — sin f + Z , 
* s a 6 a 


(54) 


where 


z = 


f na 

J rrr< 


(f.r) (R-rxT) dt 


and K s> k 6 are constants of integration. From Equations 48 and 50, and considering 


B 


f ,. 3 V 1 - e 2 r . . 

cos r dt e nt + sin f + 

2 2 a 


2 V 1 - e 2 a 2 


— — sin f , 


I 


n — sin f dt = 


$(■ 


cos f + — e cos 2 f 
2 


) 


we obtain 


f (2s + 3w) dt = f(2 A + 3B) dt + ~ (K 3 - e nt + 1 (2 sin f + -i- e sin 2 f\ 

J J n y 1-e 2 a 2 \ 2 / 

K 2 r 2 (l n , 1 nA v 

2 ( — e + 2 cos f + — e cos 2 f 1 + K A , 


n 


(55) 


(56) 


where K 4 is the additive constant of integration. 

The forms in Equations 30 to 39 of the disturbing terms in the right sides of the variational 
equations require the decomposition of r* , r^p, ?kpq , ••• in the moving frame r., r. x r , r ( 
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rather than in the frame r. , v. , Setting 

X = g T° + T} Rxr 0 + £ R 


(57) 


and taking 


—o -* dr na e sin f 
r *v = - 


dt VjTe* 


v*R x r° = 


a 2 


into consideration, we deduce from Equation 41 that 


f = (s + 2w) r - 


na e sin f 


J(2 s + 3w) dt , 


na 2 i 1 - e 2 f 

V J(2s + 3w ) < ^ t * 


(58) 

(59) 


The decomposition in Equation 57 was suggested by Popovic (References 4 and 5) and inde- 
pendently by the author (Reference 6). In Popovic’s work, the perturbations of the first and of 
the second order in ??, £ are determined. 

The disturbing vector F also can be decomposed along r , Rxr , and R. Perhaps from the 
computational standpoint this decomposition is the simplest one; it appears in several theories 
of the general planetary perturbations, either directly or as an intermediary step. Let us put 

S = Fxr , 

T = F.Rxr . 


From 


n dr -> na 2 /l-e 2 


V r d"£ 


Rx r , 


we deduce 


— _ - na 2 y 1 - e 2 n dr 

F . vxR= — — S-^T. 


Substituting this value into Equation 53, we obtain 

*■/ 


na 2 /IT e‘ _ n dr 


s - 7^ t ' +2vT 


rxR 


dt ; 
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and, replacing v by its decomposition as given above, we have-after some easy vectorial 
transformations- 


A = 



+ NT) dt, 


( 53 ’) 


where 


M = 


na 


fl^ 



r r 

sm 

a a 


(f-f) . 



r 


na 

(l - e 2 ) 3/2 




1L sin {f-i) 


2 a (l - e 2 ) r r 
r a "a 


cos 



M is a sine series in l and l, and N is a cosine series in the same arguments. To obtain these 
series, we have to obtain the development of 

— , — , — cos f , and — sin f. 
a r a a 


Equation 49 can be written as 


B = 


na T 

/u? 


dt . 


(490 


The computation of integrands in Equations 49 T , 53 T , and 55 requires also the development of 
some other expressions: for example, the odd powers of 1/p. A ; the powers of a k /r k ; the scalar 
products r. . 7. and T. •R i ; and the triple products R. ■ r.xT. . It seems that the simplest 
and easiest way to obtain all these developments is by means of the single and double harmonic 
analyses. The formulas 


7 . - A . ( cos e.-e.) +B. sine., 

i i \ 11/ i i 

e . - e . sin e . = % . , 

ii ii 

A. = a. P. , B. = a. x/i-e 2 Q- 

1 1 1 1 if 1^1 

serve as a start. The representation of p 2 . in terms of the vectorial elements is very useful 
in performing the double harmonic analysis. 


DETERMINATION OF THE CONSTANTS OF INTEGRATION 


Each perturbation 


rkp 


fTkpq 
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introduces six constants of integration. We designate them by 


K k ., K k P, K k ? q , 

J 1* ]> J 1 


(i, k, p, q, - 1, 2, •**, n; j - 1, 2, 6) 


correspondingly. Consequently, the problem of determining the general perturbations by develop- 
ing them into power series with respect to the masses introduces an infinite set of constants of 
integration. Of course, with the increasing number of the upper indices these constants decrease 
rapidly. 

In the final expressions for the perturbations, the constants are combined to form a set of 
only 6n independent constants; but this is not seen explicitly in a numerical planetary theory. At 
each step the constants must be determined separately from some additional conditions imposed 
on the elements or from the initial position and velocity vectors. Two types of elements are 
being commonly used: the mean elements, and the elements osculating at the initial moment of 
time. In the case of the mean elements the perturbations of the true longitude in the orbital plane, 
as defined by these elements, shall not contain the constant term, the purely secular term, and 
the terms with periods equal to the period of revolution of the planet. The perturbations in the 
"third coordinate," normal to the orbital plane, shall not contain the terms with periods equal to 
the revolution of the planet. The values of the mean elements, however, are not unique. They 
depend on the choice of the eccentric, true, or mean anomaly as the basic independent variable 
to be used in developing the perturbations of a given planet. 

Hansen (Reference 7) made use of the eccentric anomaly in his theory of minor planets. Hill 
(Reference 2) preferred the true anomaly. In both choices we gain speed of convergence of series 
giving the perturbations of the first order, in comparison with the choice of the mean anomaly. 
However, the road to computation of the higher order perturbations of the whole planetary system 
by either of these two choices will be blocked so effectively that all the gains in the first approxi- 
mation appear to be negated by the difficulties encountered in computing these higher order per- 
turbations. For this reason the use of the undisturbed mean anomalies t 19 <i 2 » ' ' ' » and 
consequently of the universal variable, time, is highly recommended in the planetary theories. 

This has been done already by Hansen and by Hill in their theories of Jupiter and Saturn. In con- 
nection with this statement we say that the elements are mean if there are no terms of the form 

K 0 , Kjt, K< c > cos K< s > sin l 

in the perturbations of the true longitudes with respect to the undisturbed orbit planes and if there 
are no terms of the form 


K< c > cos K< s > sin l 


in the "third coordinates" £. 

We express the perturbations in the true longitude in terms of the perturbations along 7° , 
r x 7°, and R . Taking into account the equations 
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B r -n B r o -o 

- r u , — - = r Rx r , 


Hi = o, _jLL= 5x?°, HI 

3r 2 3r BX BX 2 


o, 33 - = o, - 3 — r 

Br 2 BX Br BX 2 


= - r Rxr c 


we can write 


fsr — — r SX 2 

- - Sr 8X 2 + ••• 

7° + 

(r 8X + Sr SX r 8 X 2 + * * 

] Rx?° + SR 

V 2 

2 

1 

\ 6 

1 


f = Sr - — r SX 2 - i-Sr S\ 2 + •• 

■> o o 


T 7 = r 8 A + Sr S A. - — r S \ 3 + .... 


Solving these two last equations with respect to Sr, s\ gives 


£ C T l 2 

Sr = <f + 

2r 2 r 2 


sx =ii-i2Z-lH. + HjL + 

r r 2 3 r 3 r 3 


Substituting 


i = (f k +# p + £ q ) + (^ kp + f pq + f qk ) 


= (^? k + tjP + 17* 1 ) + (r) k ’ > + r]P'< + 77’ k ) + + 


SX = X k + X kp + X k P < ’ 


into Equation 63 gives: 


\ k = ZZ-, 


T? k P _ < fP 7? k + 77P 

r r 2 


L 



( 66 ) 


Putting 


A.kpq = _ J_ ftk ^pq + £p r ( qk + ftp} + — (ft ft ft + ft ft ft + ft ft ft) 

r r 2 r 3 

_ _L (<ftP v q + ^pq ft + ftk ft) _ JL ft ^p ^q _ 


a? n. 1 / 1 — e ? f » 

n, fT^lf W k =- 1 ^ j (2A k + 3B k ) dt , 


(67) 


a 2 n. t/ 1 - e 2 f 

n i ]/l - e* WJp = i — l (2 A kp + 3B kp jdt - 


f? ^ ft 


(68) 


n i V 1 - e T W kpq = - 3i ^ — 1" (2 A kpq + 3B kpq ) dt y (# k ft q + ft k + ,f q ftp) 

r i J r i 

+ -jr- (ft ft ft + ft ft ft + ft ft ft)--L(ft p ft + ft q ft + ft k ft)— ft Vi Vi, m 

substituting these values into Equations 64 to 66, and taking 59 into consideration, we obtain 




n . y 1 - e? 


ft ft) (- «5, * «■ * - 7 - 7 =? ( 2 •" •*" ft 


K k ; V 1 - e? 


tl-ef / 1 , \ a 

— ( 2 cos f . + — e . + — e . cos 2 f . ) + — - Rk . + W k 

n. \ 1 2 1 2 1 V r 2 -*i x 


(70) 


\kP 


fefftft) (- « kr) • 


Kft ( , 

- (2 sin f. + — e. sin 2f 

n. YUe* V 121 


■) 


K kp Vl-e? 


i ~JT (2 cos f. + ie. +|e. cos 2f.j +-j ft p + W kp . 


(71) 


A.kp q 


n. yifte 7 


1 \ 1 


if -r (- K $r + e i K i? q ) n i t 


K k Pq 
1 i 


( 2 


2 sin f. + — e. sin 2f 

n. / 1 -e 2 ' 1 ^ 


0 


Rkp q V 1-e 2 / 1 , . ... 

21 - (2 cos f, +1 e. + ±e. cos 2 f . ) +-} K kpq +W kpq . 


\ a? 

fi ) + 7? 


( 72 ) 
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We shall make use of the developments 



cos qf. 


sin qf. 


\ eg? + CP^ cos + <=§< cos 24. + 


SP? sin 4. + SP? sin 24. + 


The coefficients in these developments are computed either by some analytical classic procedure 
or by means of a harmonic analysis if the eccentricity is not too small. 

The terms of the form 


K 0 ., K. t, K< c > cos 4, K< s > sin 4 (73) 

must be absent in the developments of Equations 70 and 71. We have, keeping only the substantial 
terms, 


w i = a oi + ^oi ", t + a ii c ° s + p\i sin + ••• . 

W^ p = a^P + /3 ^p n . t + aJP cos 4 . + /S *p sin + • • • , 

W kpq = a kpq + flkpq n . t + a!?? 1 ! COS 4 . + /3>?P t l sin 4 . + ... 

1 o l ^ 0 1 l li i'll i 

In the process of computing the coefficients a and J3, the machine rejects automatically all the 
useless terms unless we decide to obtain a complete development of the perturbations in the true 
longitude. The conditions for the absence of the terms 73 in the right sides of Equations 70 and 
7 1 lead to the equations 


(74) 

(75) 

(76) 


+ 


+ 


K 2 Vl-e 2 
n 

n Vl-e 2 
K 2 /l - e 2 


(<r +i e +i c o 2 ) + i k 4 c o 2 '° + “o 

(2SJ' 1 +|e S?' 2 ) +/8 1 =0. 

(acj ' 1 tie cj’ 2 ) + C T 2 '°K 4 + <*j = 0 


= o, 


(77) 

(78) 

(79) 


+ 2^ C °' 2 ' 0 (- K 3 + eK i ) + ^o = 0 - 

Separating in Z\, Z*p, z^po, • • • the terms with the argument 4, we have 


(80) 
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Z J = c k £ cos ^ + s k . s in / C + ••• t 
Z kp = c kp cos + s kp sin '£+•••, 

iii 1 1 ’ 


Z kpq = c k ? q cos f. + s k ? q sin f- + ■*• (i, k, p, q, *•* - 1, 2, n) ; 

and the conditions for the absence of terms with the argument l in £ k , £ kp , £ kpq , * * * lead to 

K s Cl 1,1 + Cj = 0, (81) 

K 6 sp + s, = 0. (82) 

The lower index i and the upper indices k, p, q, • • • are omitted in Equations 77 to 82. Evi- 
dently, this is not causing any ambiguity. 

CONCLUSION 

The results given in this article represent the extension and completion of the results given 
in the author's previous articles on this subject. The theory given here also can be considered 
as a modification and generalization of Hill's planetary theory, with the latter's inconveniences 
removed. The interdependent constants of integration peculiar to Hill's theory do not appear in 
the present exposition. The solution is given in a form which permits us to write immediately 
the differential equation for the general perturbations proportional to any prescribed product of 
masses. Moreover, the vectorial formalism permits penetration into the structure of higher 
order effects without great difficulty. Programming also is facilitated by the repetition of the 
homogeneous operations. 

The formula of Equation 57 permits us to obtain easily the decomposition of r?, r^ 3 , • • 
along the axes of the inertial systems if it is considered necessary. 

On the basis of experience obtained at Goddard Space Flight Center, computing the pertur- 
bations of a given order for one planet might be expected to require only a few minutes. The 
simplicity of the methods for the general perturbations in the position vectors suggests that such 
methods will constitute one of the principal approaches to the problem in the not too distant future. 


(Manuscript received November 1, 1964) 
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Appendix A 

Basic Notations 


D. 

1 

D. . 

j i 

f 

m . 

i 

R. 

i 

r . 
r . 

i 

Sr. 

1 

7 a 

i 

r a .P 

i 

i 

V. 

1 


/*? 


V. exp (S r. *v.) 

V. exp * V. ) = -V. exp (Sp. * V.) 

gravitational constant 

mass of the i th planet; the mass of the sun is put equal to one 
unit vector normal to undisturbed orbit plane of the i th planet 

l^i I 

undisturbed position vector of the i th planet 
perturbations in position vector of the i th planet 
perturbations in r i proportional to m a 
perturbations in r . proportional to m a 
perturbations in r . proportional to m a m y 
del operator with respect to r. 
dir . 

i 

dt 

f (i + m i) 


p ki 

= IpJ 

^ki 

= 

s ?ki 

= S r -S r. 

k i 

Pli 

= 

k l 

7>l\ 

= rlP-7 a/3 

k i 

?kf 

= 7 a/3y - r a ' /Sy 

k i 


NASA-Langley, 1965 G-606 


21 




"The aeronautical and space activities of the United States shall be 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Adtninistration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, Uata compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


a t- 



